Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  809 1697 (N31E35 L=229)  1203    196        3 ( 730/1650- 875/1740)   0.0 E42-E28 N27-N33
      2  641 1650 (N23E44 L=221)  1436    228        4 ( 563/1553- 713/1724)  90.0 E46-E40 N15-N30
      3 1312 1415 (N26W07 L=272)  2167    505        1 (1176/1195-1372/1589) -51.0 E04-W14 N13-N34
      4 1191 1217 (N11W05 L=269)   322    157        0 (1119/1166-1245/1258) -14.0 W00-W08 N09-N12
      5 1628 1042 (N13W36 L=300)   399    115        0 (1612/ 994-1640/1101) -51.3 W33-W37 N10-N15
      6  274 1039 (S20E51 L=213)   858    183        1 ( 218/ 959- 338/1084)  45.0 E56-E45 S26-S16
      7  841  997 (S10E10 L=255)   266    131        0 ( 819/ 946- 885/1030)  80.5 E11-E07 S13-S07
      8  653  920 (S19E19 L=245)   451    204        0 ( 620/ 861- 702/ 978) -81.9 E21-E15 S21-S14
      9  286  916 (S27E49 L=216)   448     87        1 ( 243/ 903- 338/ 932)  11.3 E53-E43 S29-S25
     10 1763  792 (N03W54 L=319)   689    145        1 (1706/ 743-1805/ 851)  63.4 W50-W57 S01-N07
     11  610  455 (S50E14 L=250)  1704    250        3 ( 507/ 384- 707/ 542) -15.9 E25-E03 S54-S46
     12 1213  486 (S30W28 L=293)   226     58        0 (1189/ 463-1243/ 501)  63.4 W26-W29 S31-S27
  

    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