Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  716 1723 (N41E31 L=297)   838     93        1 ( 665/1675- 758/1756) -24.4 E36-E25 N36-N43
      2  592 1714 (N40E42 L=285)   421     40        0 ( 574/1685- 609/1740) -63.4 E44-E39 N37-N41
      3  436 1574 (N28E49 L=278)   588     72        0 ( 395/1555- 476/1590)  -6.3 E53-E44 N27-N29
      4 1858 1282 (N18W63 L= 30)   554    121        0 (1841/1241-1892/1315)  45.0 W61-W69 N15-N20
      5  297 1158 (N01E50 L=277)   462     71        0 ( 265/1125- 331/1190)  26.6 E53-E47 N00-N03
      6  467  840 (S18E36 L=291)   260     74        0 ( 457/ 792- 483/ 876)  68.2 E37-E34 S20-S15
      7  293  751 (S23E54 L=274)   752    129        1 ( 252/ 694- 344/ 792)  21.0 E60-E47 S26-S20
      8 1737  498 (S31W67 L= 34)  2431    212        4 (1663/ 446-1809/ 577)  39.3 W59-W75 S34-S24
      9  610  508 (S39E30 L=297)  1221    197        1 ( 505/ 471- 689/ 544)  -7.1 E38-E22 S41-S36
     10 1228  423 (S42W21 L=348)   426    103        1 (1208/ 395-1267/ 472) -33.7 W18-W24 S43-S38
     11 1591  425 (S38W55 L= 22)   473     63        0 (1576/ 390-1609/ 458) -51.3 W52-W57 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