Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1917 (N64W45 L= 37)   467     48        0 ( 950/1890- 986/1937)  24.8 W39-W52 N60-N66
      2  813 1839 (N64W17 L=  9)  2337    361        3 ( 686/1743- 931/1913)   7.8 E03-W41 N58-N67
      3  912 1616 (N42W10 L=  2)  1458    335        1 ( 763/1589-1073/1649) -20.9 E01-W20 N35-N45
      4 1681 1543 (N15W60 L= 52)   526     64        1 (1642/1522-1713/1566)   0.0 W55-W65 N14-N16
      5 1196 1500 (N26W23 L= 15)   517     86        1 (1140/1462-1228/1536)   0.0 W19-W26 N24-N28
      6 1828 1429 (N05W69 L= 61)  1602    106        2 (1799/1374-1855/1477)  50.2 W65-W75 N01-N07
      7  376 1417 (N43E39 L=313)   456     44        0 ( 347/1393- 410/1435)   0.0 E42-E35 N41-N44
      8 1500 1307 (N08W35 L= 27)   704    168        0 (1462/1238-1536/1384)  83.7 W33-W36 N03-N12
      9 1784 1227 (S04W54 L= 46)   310     22        1 (1766/1208-1805/1246) -33.7 W52-W56 S05-S02
     10 1170 1088 (N04W10 L=  1)   489    125        0 (1129/1032-1215/1143) -81.9 W08-W11 N00-N07
     11  162 1024 (N22E63 L=288)  4379    272       10 ( 100/ 933- 243/1142) -49.8 E74-E54 N16-N29
     12  397  911 (N12E41 L=310)   960    211        2 ( 324/ 855- 474/ 987) -48.0 E45-E36 N06-N17
     13 1170  790 (S12W03 L=354)   523    147        0 (1141/ 724-1199/ 871)  54.5 W00-W05 S15-S08
     14  153  712 (N04E75 L=277)  3607    195        5 ( 114/ 640- 197/ 801) -56.3 E83-E69 N00-N09
     15 1743  588 (S42W43 L= 35)   700    109        1 (1691/ 568-1790/ 611) -18.4 W36-W51 S45-S38
     16  552  512 (S14E43 L=309)   806    167        0 ( 525/ 438- 587/ 593)  36.0 E48-E37 S18-S10
     17 1571  510 (S42W25 L= 16)   484     81        0 (1536/ 480-1615/ 542)  14.0 W21-W29 S43-S40
  

    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