Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  702 1339 (N30E17 L=  6)  1497    385        3 ( 541/1286- 887/1389)  -2.2 E29-E03 N27-N31
      2 1318 1218 (N13W21 L= 44)   647    160        2 (1249/1164-1387/1250) -45.0 W17-W24 N08-N15
      3  177 1044 (N16E64 L=319)   892    193        0 ( 139/1011- 222/1090) -21.8 E68-E58 N13-N18
      4 1670  975 (S07W41 L= 65)   954    234        0 (1625/ 865-1712/1083)  53.1 W37-W46 S13-S01
      5  684  730 (S07E25 L=358)  1606    350        5 ( 603/ 582- 729/ 882)  57.1 E34-E21 S15-N02
      6  943  796 (S06E08 L= 15)   639    198        2 ( 856/ 762-1038/ 826) -20.0 E13-E02 S09-S04
      7 1758  749 (S24W48 L= 71)  1029    129        3 (1705/ 706-1800/ 783)  -6.3 W43-W52 S26-S21
      8 1565  715 (S22W31 L= 54)   353    101        0 (1536/ 674-1582/ 757) -71.6 W29-W32 S24-S18
      9  853  482 (S25E20 L=  3)   643    247        3 ( 760/ 403- 914/ 553)  23.2 E28-E14 S30-S21
     10  652  306 (S37E43 L=340)  1367    205        2 ( 611/ 256- 691/ 376)  68.2 E48-E40 S41-S30
  

    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