Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  946 1583 (N25E21 L=334)   513    248        0 ( 924/1481- 978/1690) -60.3 E24-E16 N19-N33
      2 1181 1535 (N29E05 L=350)   387    141        0 (1133/1473-1223/1565) -20.6 E08-W00 N26-N30
      3  511 1386 (N03E41 L=315)  1609    360        2 ( 360/1308- 611/1504)  55.5 E49-E36 S04-N12
      4 1267 1286 (N16W07 L=  2)   558    216        2 (1207/1197-1328/1356) -28.6 W01-W12 N12-N19
      5  216 1199 (S13E61 L=294)   682    160        2 ( 189/1133- 242/1280)  90.0 E64-E58 S17-S08
      6  266 1166 (S14E54 L=301)   323     41        1 ( 246/1146- 290/1185)   0.0 E56-E51 S15-S13
      7 1052 1124 (N02E01 L=354)   580    207        0 ( 977/1088-1150/1167)  11.3 E06-W04 N00-N03
      8  761 1015 (S12E15 L=341)   185     51        0 ( 738/ 993- 788/1040)  63.4 E15-E13 S13-S09
      9 1514  849 (S01W33 L= 28)   213     72        1 (1493/ 821-1534/ 881) -71.6 W31-W34 S02-N01
     10  899  690 (S27W02 L=357)   207    104        0 ( 868/ 663- 943/ 725) -11.3 W00-W05 S27-S25
     11  610  568 (S43E15 L=340)  1143    212        3 ( 508/ 536- 713/ 603)  21.3 E24-E06 S47-S38
     12 1284  300 (S38W47 L= 42)   751    178        0 (1233/ 276-1363/ 347)  15.3 W39-W51 S40-S34
  

    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