Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1562 1736 (N30W66 L= 74)  1414    159        0 (1522/1686-1629/1769)   0.0 W58-W83 N27-N31
      2 1828 1179 (S08W57 L= 65)  1848    274        7 (1784/1052-1862/1287)  61.7 W53-W62 S14-N00
      3  450 1226 (N27E32 L=336)   274     89        0 ( 429/1189- 486/1267) -56.3 E34-E29 N24-N30
      4  219 1146 (N27E55 L=313)   378     64        0 ( 186/1129- 253/1165)   0.0 E59-E51 N26-N27
      5 1344 1029 (S04W18 L= 26)  1760    571        5 (1157/ 849-1469/1162) -53.6 W08-W22 S16-N05
      6 1612 1096 (S07W37 L= 45)   566    122        2 (1570/1032-1653/1146) -82.9 W34-W39 S10-S02
      7 1697 1045 (S11W42 L= 50)  1125    233        5 (1643/ 956-1737/1135)  38.7 W36-W46 S15-S06
      8  410 1114 (N21E37 L=331)   303    131        1 ( 346/1098- 472/1129) -10.3 E42-E31 N19-N22
      9  766 1051 (N10E14 L=354)   466    161        1 ( 689/1031- 843/1078) -29.1 E18-E09 N08-N13
     10 1762  813 (S27W45 L= 53)  1415    249        3 (1678/ 781-1874/ 861)  -8.1 W37-W58 S30-S22
     11 1572  731 (S27W28 L= 36)   314     87        0 (1561/ 687-1587/ 778)  71.6 W26-W28 S29-S23
     12  917  723 (S11E13 L=355)   546    141        3 ( 856/ 679- 965/ 767)   7.1 E17-E09 S13-S09
     13 1590  573 (S38W28 L= 36)  2236    299        7 (1520/ 445-1663/ 697)  53.7 W22-W35 S46-S31
     14  979  586 (S21E13 L=355)   796    241        1 ( 893/ 514-1035/ 692)  30.3 E19-E07 S24-S15
     15  511  593 (S10E42 L=326)   440    131        0 ( 435/ 573- 563/ 623)   0.0 E48-E37 S11-S09
     16  695  453 (S23E35 L=332)  1499    344        1 ( 556/ 363- 815/ 545)  10.9 E50-E24 S26-S19
     17  545  307 (S29E57 L=310)   765     86        1 ( 518/ 253- 580/ 362)  36.0 E62-E50 S33-S25
  

    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