Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1149 1742 (N54W14 L=164)   475     83        0 (1115/1711-1182/1767)  14.0 W10-W18 N51-N56
      2  957 1553 (N38E05 L=145)   677    200        2 ( 866/1524-1038/1585)  -9.5 E11-W01 N36-N40
      3  313 1454 (N30E61 L= 89)  1713    255        1 ( 268/1359- 371/1564) -20.6 E82-E49 N24-N36
      4 1170 1499 (N34W11 L=161)   263     67        0 (1134/1478-1212/1519)  18.4 W08-W14 N32-N35
      5  711 1252 (N18E21 L=130)   807    298        2 ( 620/1180- 800/1300)  14.0 E26-E14 N13-N21
      6  306 1230 (N16E53 L= 97)   433    112        1 ( 294/1177- 324/1286) -81.9 E54-E50 N12-N19
      7  362 1019 (N03E45 L=105)  1034    232        0 ( 302/ 936- 445/1099)  39.8 E50-E38 S02-N08
      8  316  891 (S05E50 L=100)   610     80        1 ( 284/ 866- 356/ 923)  15.9 E53-E46 S06-S03
      9  303  677 (S19E56 L= 94)  1528    190        3 ( 276/ 592- 349/ 735)  45.0 E61-E49 S25-S15
     10 1591  691 (S18W39 L=189)   392    138        1 (1518/ 677-1647/ 718)   5.2 W33-W44 S19-S16
     11  716  572 (S25E22 L=128)  1458    457        1 ( 565/ 427- 876/ 715)  34.6 E38-E09 S36-S15
     12 1002  482 (S32E02 L=148)  3169    654        7 ( 763/ 347-1227/ 690)  35.2 E22-W12 S42-S17
  

    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