Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  818 1890 (N68E02 L=327)   923     70        1 ( 798/1847- 842/1922)  61.4 E04-W04 N62-N74
      2  657 1462 (N32E19 L=310)   420     94        0 ( 609/1447- 705/1480)   0.0 E22-E15 N30-N33
      3 1163 1365 (N18W13 L=342)   971    333        2 (1035/1239-1239/1456) -63.4 W06-W16 N09-N25
      4 1804 1415 (N13W65 L= 33)  1337    152        2 (1730/1365-1863/1448)  -6.0 W56-W75 N09-N14
      5  663 1416 (N29E18 L=310)   441     84        1 ( 616/1395- 712/1433)  -8.1 E21-E14 N27-N30
      6  302 1243 (N23E48 L=281)   546    122        1 ( 283/1181- 323/1300) -76.0 E49-E46 N18-N26
      7 1170 1083 (N02W09 L=338)   303     57        1 (1145/1054-1194/1118)  14.0 W07-W11 N00-N03
      8  985  826 (S11E05 L=324)   194     84        0 ( 951/ 808-1024/ 840) -26.6 E07-E02 S12-S10
      9 1531  559 (S36W30 L=359)  2838    471        2 (1340/ 435-1709/ 686)  17.4 W14-W46 S42-S29
     10  465  497 (S24E49 L=280)  3518    523        6 ( 371/ 293- 600/ 678)  76.0 E54-E43 S39-S11
     11 1120  175 (S62E12 L=316)   434     33        0 (1101/ 158-1144/ 199) -71.6 E15-E09 S64-S58
  

    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