Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  792 1871 (N47E42 L=197)   878    173        0 ( 772/1833- 813/1911) -41.6 E46-E36 N43-N51
      2 1057 1536 (N25E07 L=232)   751    174        1 (1012/1450-1110/1596)  79.7 E08-E04 N18-N29
      3 1701 1483 (N36W41 L=280)  1511    211        2 (1633/1443-1765/1520)   4.4 W34-W47 N33-N39
      4 1381 1432 (N24W15 L=254)   940    187        2 (1324/1334-1430/1517) -45.0 W10-W19 N18-N29
      5  297 1454 (N09E59 L=179)   605    112        2 ( 266/1411- 323/1491)  78.7 E62-E57 N06-N11
      6  391 1426 (N08E49 L=190)   573    114        0 ( 358/1371- 420/1481)  74.1 E50-E47 N04-N12
      7 1763 1292 (N24W45 L=284)   319     58        0 (1738/1277-1790/1308)   0.0 W43-W48 N23-N25
      8 1689 1217 (N18W39 L=278)   999    201        3 (1655/1138-1717/1303)  90.0 W36-W41 N13-N23
      9  421 1097 (S11E39 L=200)  1620    395        3 ( 343/ 927- 488/1262)  74.1 E44-E36 S21-N00
     10 1210 1206 (N07W08 L=246)   474    127        1 (1165/1161-1248/1248) -26.6 W04-W10 N04-N09
     11  137 1173 (S08E70 L=169)   915    169        1 ( 120/1131- 162/1230) -50.2 E74-E65 S10-S04
     12  876  997 (S11E08 L=231)   183     64        1 ( 844/ 977- 920/1012)  26.6 E10-E05 S11-S09
     13  188  822 (S29E59 L=180)  1502    420        0 ( 128/ 723- 248/ 917)  40.6 E66-E52 S35-S23
     14 1371  790 (S14W25 L=264)  2842    710        3 (1027/ 725-1663/ 834)  18.0 W04-W44 S19-S06
     15  822  670 (S31E07 L=232)   616    211        1 ( 720/ 656- 926/ 689)  25.7 E13-W00 S32-S28
     16  543  599 (S40E28 L=211)  1422    373        5 ( 422/ 537- 685/ 659)  22.2 E38-E16 S45-S33
     17 1732  485 (S22W65 L=303)  2150    253        3 (1660/ 438-1801/ 533)  38.7 W58-W70 S25-S17
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y