Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  406 1590 (N58E47 L=326)  1092    187        2 ( 316/1561- 479/1626)   0.0 E66-E35 N55-N60
      2 1570 1416 (N15W45 L= 57)   565    143        0 (1541/1369-1604/1483)  60.3 W42-W50 N11-N18
      3 1155 1292 (N19W14 L= 26)  1674    560        0 ( 892/1244-1371/1339) -21.4 E01-W27 N13-N25
      4  615 1314 (N33E20 L=352)   489    146        1 ( 546/1286- 685/1334) -26.6 E25-E15 N29-N35
      5  164 1209 (N34E69 L=303)  1715    104        3 ( 111/1169- 209/1240)  -9.5 E85-E60 N31-N36
      6  859  989 (N09E10 L=  2)   854    260        2 ( 789/ 888- 915/1090)  85.2 E12-E07 N03-N15
      7 1801  953 (S17W50 L= 62)   961    164        2 (1742/ 919-1872/ 992) -15.9 W44-W58 S19-S13
      8 1291  771 (S13W10 L= 22)   440    146        1 (1238/ 708-1321/ 807) -74.1 W07-W11 S17-S10
      9  241  583 (S05E71 L=301)   566     66        0 ( 222/ 562- 267/ 614)  45.0 E74-E66 S06-S02
     10  752  550 (S15E28 L=344)   858    136        1 ( 718/ 482- 793/ 608)  41.2 E31-E23 S18-S11
     11 1480  398 (S45W17 L= 29)  1145    255        0 (1371/ 346-1591/ 448)  -2.6 W06-W28 S47-S43
     12  604  271 (S31E56 L=316)  1566    137        1 ( 563/ 217- 638/ 326)  20.0 E70-E48 S35-S27
     13 1188  243 (S49E13 L=359)   395     64        0 (1151/ 222-1220/ 263)   0.0 E16-E09 S50-S47
  

    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