Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  522 1468 (N29E32 L=171)   495    100        0 ( 488/1423- 553/1513) -45.0 E35-E29 N26-N32
      2  164 1358 (N26E69 L=134)   752     55        2 ( 135/1338- 195/1384) -15.9 E78-E63 N23-N27
      3 1715 1294 (N10W49 L=252)   315     92        0 (1694/1244-1728/1341)  80.5 W47-W49 N07-N13
      4  200 1125 (N10E58 L=145)   722     87        2 ( 167/1086- 220/1175)  78.7 E61-E55 N08-N13
      5  245 1130 (N10E53 L=150)   301     35        0 ( 229/1107- 262/1156)  45.0 E54-E51 N08-N11
      6  526 1029 (N02E31 L=172)   221     86        0 ( 485/1008- 566/1044)   0.0 E33-E27 N01-N02
      7 1137  972 (S06W06 L=209)   393     96        0 (1117/ 923-1160/1024) -80.5 W05-W07 S08-S02
      8 1314  952 (S08W17 L=220)   212     89        0 (1303/ 905-1327/ 999)  78.7 W16-W17 S10-S05
      9  503  722 (S16E36 L=167)  1619    414        6 ( 379/ 560- 586/ 862)  36.0 E50-E28 S24-S08
     10  431  773 (S12E41 L=162)   694      0        0 ( 359/ 734- 490/ 804) -16.7 E46-E36 S15-S10
     11  356  714 (S15E48 L=155)   516     49        1 ( 337/ 675- 374/ 760)  51.3 E50-E46 S17-S12
     12  180  662 (S17E71 L=132)  1649    132        5 ( 143/ 617- 212/ 714)  28.6 E80-E65 S19-S13
     13 1591  577 (S33W39 L=242)   737    191        0 (1501/ 531-1679/ 614)  15.9 W31-W47 S35-S30
     14 1147  535 (S33W05 L=208)   227     62        0 (1107/ 521-1178/ 548)   0.0 W01-W06 S33-S31
     15  534  464 (S33E42 L=161)   479     52        0 ( 501/ 442- 564/ 488)  15.9 E45-E38 S34-S31
     16 1241  454 (S39W12 L=215)   737    223        1 (1136/ 408-1344/ 478)   7.1 W03-W19 S42-S37
     17  666  289 (S47E41 L=162)   517     81        0 ( 630/ 256- 694/ 326)  22.6 E48-E36 S49-S44
  

    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