Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  927 1801 (N45E35 L=348)   338     79        0 ( 904/1767- 958/1830) -74.1 E36-E32 N41-N48
      2  535 1672 (N25E55 L=328)   965    103        3 ( 489/1635- 579/1697)   0.0 E61-E48 N23-N26
      3  387 1495 (N11E56 L=326)  1058    177        3 ( 338/1433- 424/1557) -45.0 E61-E50 N07-N15
      4 1256 1395 (N24W05 L= 28)   286     73        0 (1218/1369-1297/1414)  -9.5 W02-W08 N23-N25
      5 1399 1115 (N11W20 L= 43)   281     89        0 (1376/1065-1413/1161) -59.0 W17-W21 N08-N13
      6 1835  962 (N14W59 L= 81)   887    150        0 (1789/ 891-1862/1005) -21.8 W52-W62 N09-N15
      7  308  825 (S31E47 L=335)   500     99        1 ( 273/ 783- 337/ 863)  40.6 E51-E44 S34-S28
      8 1091  611 (S25W15 L= 37)   218     52        0 (1062/ 594-1119/ 628) -14.0 W12-W16 S26-S24
      9 1111  455 (S35W22 L= 45)   759    154        2 (1036/ 418-1172/ 501)  -9.5 W15-W27 S35-S32
     10 1564  442 (S23W56 L= 79)   860    142        0 (1505/ 418-1632/ 475)  26.6 W50-W61 S24-S19
     11 1220  261 (S45W44 L= 67)   421     74        0 (1200/ 245-1260/ 290)   7.1 W40-W48 S46-S42
  

    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