Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1580 1526 (N33W46 L= 40)   560    119        0 (1532/1499-1632/1558) -39.8 W42-W50 N30-N35
      2  736 1337 (N25E19 L=335)  3745   1127       11 ( 330/1095-1034/1552) -27.8 E57-W00 N08-N39
      3 1032 1465 (N32W02 L=356)   212     59        0 (1004/1447-1068/1477)  11.3 W00-W05 N31-N33
      4 1444 1415 (N27W32 L= 25)   713    205        4 (1363/1357-1536/1439) -26.6 W25-W37 N22-N28
      5 1674 1225 (N13W46 L= 40)   905    229        0 (1590/1142-1760/1294) -39.8 W40-W53 N07-N18
      6  324 1176 (N15E49 L=304)   902    207        2 ( 302/1091- 348/1285) -76.0 E51-E46 N09-N21
      7  155  927 (S01E70 L=284)   788    177        2 ( 121/ 906- 194/ 944)   8.1 E77-E63 S03-N00
      8 1391  859 (S08W23 L= 16)   458    120        0 (1332/ 818-1444/ 902) -32.0 W18-W26 S10-S04
      9  894  693 (S16E09 L=344)  1125    365        3 ( 753/ 557- 996/ 832)  43.4 E20-E02 S24-S07
     10 1542  759 (S15W34 L= 27)   247     55        0 (1512/ 745-1563/ 778)   0.0 W31-W35 S15-S13
     11 1077  501 (S30W02 L=355)  1803    441        5 ( 932/ 371-1254/ 602)  34.8 E10-W13 S39-S23
  

    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