Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1570 1702 (N27W63 L=139)   595     62        1 (1544/1674-1594/1727)   6.3 W59-W68 N25-N27
      2 1394 1608 (N27W41 L=117)   690    118        1 (1338/1571-1450/1642) -60.3 W37-W43 N23-N30
      3 1708 1539 (N14W63 L=138)  2271    215        6 (1642/1479-1769/1617)  10.8 W54-W75 N09-N17
      4 1344 1412 (N17W29 L=105)  1250    401        0 (1191/1369-1532/1505) -45.0 W22-W40 N10-N26
      5  981 1384 (N25W07 L= 82)  2283    639        3 ( 719/1251-1220/1444) -30.7 E10-W18 N11-N32
      6  856 1241 (N20E04 L= 71)   455    200        1 ( 757/1227- 950/1256) -31.0 E10-W00 N16-N22
      7 1588 1129 (S04W36 L=111)   442     80        1 (1567/1080-1604/1177)  53.1 W33-W37 S07-S01
      8 1876 1008 (S19W59 L=134)  1587    160        1 (1812/ 961-1946/1043)  -2.9 W50-W70 S20-S16
      9  787  934 (N04E15 L= 60)   347    116        1 ( 753/ 887- 831/ 977)  26.6 E17-E11 N02-N05
     10 1200  870 (S09W06 L= 82)   348    173        2 (1155/ 823-1225/ 949)  50.2 W02-W08 S11-S04
     11  806  552 (S18E25 L= 51)  3009    715       10 ( 643/ 413-1047/ 770) -55.7 E34-E12 S28-S03
     12  279  719 (N03E57 L= 19)  1439    165        2 ( 207/ 673- 351/ 760) -15.9 E65-E51 S01-N04
     13  609  661 (S07E33 L= 42)   627    140        2 ( 577/ 604- 644/ 729)  50.2 E36-E30 S10-S03
     14 1016  626 (S19E10 L= 65)   368    105        1 ( 963/ 602-1066/ 653) -18.4 E13-E07 S19-S17
     15  288  533 (S08E66 L= 10)  4062    215       11 ( 240/ 447- 337/ 637)  47.5 E78-E59 S13-S01
     16  458  516 (S13E50 L= 26)  1182    227        3 ( 384/ 468- 560/ 591) -41.6 E53-E42 S16-S07
     17  704  298 (S32E45 L= 31)  4015    347       10 ( 638/ 131- 770/ 441)  34.4 E74-E37 S45-S21
  

    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