Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1497 1605 (N27W45 L=238)  1162    215        2 (1426/1564-1561/1639)   8.1 W37-W51 N25-N28
      2  842 1562 (N37E02 L=192)   823    129        0 ( 768/1530- 901/1589) -18.4 E07-W02 N34-N39
      3 1348 1560 (N27W32 L=226)   273     52        0 (1332/1525-1365/1584)  45.0 W30-W33 N25-N29
      4 1598 1358 (N11W43 L=237)   325     66        0 (1569/1341-1637/1376) -26.6 W40-W45 N09-N12
      5  586 1277 (N23E23 L=170)   841    201        0 ( 500/1227- 663/1332) -36.0 E29-E18 N19-N28
      6 1280 1269 (N11W19 L=213)   276     63        1 (1249/1243-1319/1296) -45.0 W17-W21 N09-N13
      7  926 1240 (N15E02 L=192)   739    232        2 ( 857/1166- 982/1291)  35.5 E07-W01 N11-N18
      8  216 1225 (N26E55 L=139)   540    110        0 ( 188/1195- 260/1245)  15.9 E57-E50 N24-N27
      9 1767 1218 (S00W53 L=247)   362     90        0 (1739/1202-1802/1235)  -9.5 W50-W56 S01-N00
     10  774  422 (S31E29 L=164)   475     97        0 ( 729/ 396- 816/ 454)  14.0 E33-E25 S32-S29
     11  960  211 (S52E28 L=166)   383     22        0 ( 936/ 194- 983/ 227) -56.3 E29-E25 S54-S50
  

    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