Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1297 1730 (N54W17 L= 83)   324     48        0 (1268/1717-1325/1744)   0.0 W14-W20 N52-N55
      2  593 1476 (N25E36 L= 30)  1439    331        2 ( 504/1392- 646/1606) -82.4 E42-E31 N19-N34
      3 1498 1387 (N29W30 L= 96)  2887    702        5 (1233/1173-1691/1529) -27.2 W10-W45 N16-N37
      4 1345 1199 (N15W19 L= 85)   390    151        1 (1314/1146-1390/1272) -54.5 W16-W22 N12-N19
      5 1531 1197 (N17W32 L= 98)   310     98        0 (1514/1145-1544/1244) -80.5 W30-W33 N13-N19
      6  432  781 (S20E39 L= 27)   681    143        1 ( 411/ 715- 463/ 842)  76.0 E40-E35 S23-S15
      7 1448  652 (S18W32 L= 98)   781    205        2 (1367/ 590-1528/ 709)  45.0 W27-W37 S22-S13
      8 1711  652 (S16W55 L=121)   475     54        0 (1693/ 623-1732/ 677)  76.0 W53-W57 S17-S13
      9 1714  569 (S21W60 L=125)  1370    147        4 (1655/ 516-1782/ 617)  45.0 W53-W68 S24-S17
  

    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