Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1043 1910 (N53E43 L=293)   563     64        1 (1018/1880-1070/1936) -42.0 E49-E39 N48-N57
      2  612 1659 (N19E47 L=289)  2799    473        8 ( 546/1452- 693/1861) -48.0 E66-E38 N05-N35
      3 1095 1580 (N28E13 L=323)   237     66        0 (1075/1550-1114/1612) -31.0 E15-E10 N26-N29
      4  484  773 (S35E28 L=308)  1805    357        1 ( 330/ 704- 637/ 831)  31.6 E44-E17 S43-S27
      5 1136  778 (S16W13 L=350)   416    107        1 (1084/ 750-1183/ 798)  18.4 W10-W16 S17-S15
      6  972  666 (S27W08 L=344)   722    148        1 ( 884/ 647-1051/ 689)  31.0 W02-W12 S30-S24
      7 1688  654 (S05W54 L= 30)   491     79        1 (1638/ 637-1731/ 672)  20.6 W49-W57 S07-S03
      8 1306  607 (S20W30 L=  6)   732    115        2 (1270/ 570-1347/ 648)  80.5 W26-W32 S23-S17
      9 1420  429 (S26W47 L= 24)   386    108        0 (1371/ 413-1476/ 444)  39.8 W43-W50 S28-S23
     10  536  375 (S64E20 L=317)  1634    167        1 ( 461/ 345- 618/ 403)  19.1 E34-E08 S68-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