Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  661 1830 (N40E60 L= 60)  1456    122        1 ( 629/1777- 701/1880) -10.6 E82-E50 N36-N43
      2  882 1644 (N31E26 L= 94)   563    162        0 ( 853/1578- 922/1708) -52.1 E29-E20 N27-N36
      3  330 1555 (N14E66 L= 54)  3101    196        4 ( 265/1489- 381/1651) -85.2 E73-E57 N08-N20
      4  421 1605 (N19E60 L= 60)  1014     95        2 ( 390/1572- 457/1644) -16.7 E64-E54 N16-N21
      5  233 1329 (S01E64 L= 56)   607    157        1 ( 187/1259- 255/1389)  82.9 E67-E61 S06-N02
      6  979 1086 (S00E04 L=116)  1204    343        2 ( 936/ 915-1046/1223) -51.8 E08-W03 S08-N06
      7 1783 1089 (N19W50 L=170)  1632    303        1 (1692/ 972-1854/1205) -25.3 W39-W59 N13-N23
      8  475 1118 (S09E36 L= 84)   371    118        0 ( 423/1100- 530/1156)  35.5 E39-E32 S10-S05
      9 1566  612 (S13W45 L=165)  1476    297        3 (1431/ 544-1696/ 664)  36.9 W37-W53 S20-S07
     10 1581  412 (S24W59 L=179)  2145    201        5 (1489/ 383-1672/ 442)  14.7 W50-W70 S27-S21
     11  909  378 (S46W11 L=131)   730    124        1 ( 834/ 359- 965/ 412)  10.3 W05-W16 S48-S42
  

    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