Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  395 1491 (N31E53 L=234)   429     76        1 ( 372/1466- 430/1526) -26.6 E55-E47 N29-N33
      2 1630 1451 (N30W48 L=335)  1045    147        2 (1575/1384-1688/1503) -53.1 W43-W51 N25-N33
      3 1769 1469 (N30W67 L=354)   762     76        1 (1726/1442-1802/1492) -20.6 W60-W71 N28-N32
      4  370 1024 (N02E44 L=243)   301     54        0 ( 346/ 993- 390/1051) -45.0 E46-E42 N00-N03
      5 1682  958 (S01W45 L=332)   745    109        2 (1651/ 897-1712/ 998) -59.0 W42-W47 S04-N01
      6  154  873 (S09E70 L=217)  1507    142        3 ( 131/ 822- 175/ 933)  60.3 E76-E66 S12-S05
      7 1829  858 (S08W61 L=348)   446     50        0 (1811/ 819-1844/ 882) -53.1 W58-W62 S10-S06
      8 1769  705 (S18W57 L=344)  1125    204        2 (1703/ 665-1843/ 747)  14.9 W50-W66 S19-S15
      9 1186  681 (S18W11 L=298)   323     93        0 (1168/ 633-1202/ 726) -71.6 W09-W11 S21-S15
     10  203  686 (S21E69 L=218)  1576     99        3 ( 169/ 646- 238/ 723)  16.4 E80-E62 S23-S17
     11 1787  588 (S26W66 L=353)  1711    255        2 (1718/ 508-1849/ 691)  56.3 W59-W70 S31-S19
     12 1103  512 (S30W06 L=293)  1614    349        2 ( 940/ 460-1247/ 575)  21.8 E05-W15 S33-S25
  

    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