Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  872 1772 (N48E30 L=244)   447     46        1 ( 848/1752- 904/1789)  45.0 E32-E26 N45-N49
      2  536 1548 (N25E44 L=230)   999    163        1 ( 479/1469- 593/1628)  85.2 E46-E41 N19-N31
      3  336 1454 (N16E58 L=216)   538    179        1 ( 297/1417- 375/1487) -10.3 E63-E52 N14-N17
      4  390 1335 (N10E48 L=225)   605     97        1 ( 372/1287- 411/1387) -71.6 E50-E46 N06-N12
      5 1133 1271 (N17W03 L=277)   298     94        1 (1091/1225-1163/1313) -51.3 W00-W05 N13-N18
      6  922 1012 (S02E06 L=268)   996    298        1 ( 854/ 887-1024/1140) -61.7 E10-W01 S08-N05
      7  401  968 (S12E40 L=233)   258     74        0 ( 366/ 956- 430/ 979)  18.4 E43-E37 S13-S11
      8  682  933 (S10E20 L=254)   405    161        1 ( 631/ 869- 751/ 972)  29.7 E22-E15 S14-S07
      9 1036  820 (S12W04 L=277)   216     42        2 (1008/ 801-1066/ 843)  45.0 W02-W05 S13-S10
  

    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