Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1762 (N37E30 L=156)   292     46        1 ( 951/1733- 983/1794) -39.8 E33-E27 N34-N39
      2  413 1671 (N16E68 L=119)  1767     87        1 ( 378/1640- 462/1704)  11.3 E79-E62 N13-N18
      3 1462 1607 (N43W11 L=197)   565    149        0 (1408/1546-1508/1665) -24.4 W04-W15 N39-N46
      4 1377 1569 (N37W05 L=192)   353    117        0 (1349/1511-1403/1621) -40.6 W01-W08 N33-N40
      5  661 1280 (S02E28 L=159)  1174    320        2 ( 593/1135- 732/1407) -36.3 E36-E20 S08-N04
      6 1574 1179 (N19W28 L=215)   915    193        0 (1512/1097-1631/1236) -24.4 W22-W33 N15-N20
      7 1225 1176 (N08W07 L=193)   209     53        1 (1207/1153-1248/1208) -56.3 W05-W08 N06-N09
      8  199 1152 (S19E63 L=124)   431     64        0 ( 180/1120- 219/1188)  68.2 E65-E60 S21-S16
      9  102 1101 (S23E81 L=106)  1020     76        0 (  91/1063- 115/1140)  68.2 E85-E77 S25-S20
     10 1645  866 (N04W42 L=228)   476    123        0 (1627/ 808-1667/ 923) -39.8 W39-W45 N01-N06
     11  778  714 (S30E06 L=181)   896    291        0 ( 664/ 660- 908/ 780)  49.4 E12-W00 S36-S22
     12 1715  748 (N00W52 L=238)   325     17        0 (1692/ 734-1741/ 763)   0.0 W49-W54 N00-N01
     13 1194  732 (S17W18 L=205)  1018    310        0 (1044/ 701-1335/ 761)  33.7 W10-W26 S22-S12
     14 1257  453 (S30W35 L=222)  1033    228        0 (1161/ 432-1366/ 471)  39.3 W29-W40 S34-S25
     15  488  399 (S64E28 L=159)   455     34        0 ( 468/ 378- 511/ 426)  76.0 E31-E23 S66-S61
     16  573  372 (S62E14 L=172)   457     62        0 ( 545/ 348- 602/ 401) -20.6 E18-E09 S64-S60
     17 1487  332 (S29W61 L=248)   460     51        0 (1471/ 299-1501/ 364) -18.4 W56-W65 S30-S27
  

    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