Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1046 1646 (N41E10 L=123)   889    177        2 (1017/1573-1078/1700) -45.0 E14-E05 N35-N44
      2  451 1519 (N22E50 L= 84)  1076    313        1 ( 382/1414- 513/1656) -86.8 E53-E46 N14-N32
      3 1461 1574 (N43W25 L=159)   392     91        0 (1431/1554-1499/1592)  23.2 W22-W29 N40-N44
      4  634 1506 (N24E35 L= 99)   385    122        1 ( 611/1459- 652/1567) -54.5 E36-E31 N21-N28
      5  337 1527 (N21E61 L= 72)   463     69        0 ( 316/1507- 358/1547) -18.4 E63-E57 N19-N22
      6 1351 1455 (N32W16 L=149)   300     60        1 (1332/1423-1371/1485) -78.7 W14-W17 N29-N34
      7  788 1339 (N16E20 L=114)   218     91        1 ( 736/1320- 827/1359)  31.0 E22-E16 N13-N17
      8 1011 1246 (N13E04 L=129)   286     69        0 ( 978/1228-1051/1264)  14.0 E05-E01 N11-N14
      9 1091 1091 (N05W03 L=136)   313    105        1 (1052/1033-1114/1145) -56.3 W00-W05 N01-N07
     10  649  935 (S11E22 L=111)   332    127        1 ( 584/ 919- 705/ 955)   7.1 E26-E18 S11-S08
     11 1010  845 (S11W02 L=135)   835    203        3 ( 948/ 767-1047/ 922) -66.0 E02-W05 S15-S06
     12  464  563 (S38E36 L= 98)   445    111        1 ( 423/ 538- 500/ 593)  32.0 E39-E31 S40-S35
  

    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