Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  479 1540 (N29E49 L=213)  1906    333        3 ( 437/1423- 565/1651) -36.9 E61-E37 N21-N36
      2 1453 1520 (N38W29 L=290)   738    163        1 (1396/1461-1513/1567) -33.7 W24-W33 N33-N41
      3 1163 1203 (N14W07 L=268)   243    117        0 (1106/1168-1205/1220) -15.9 W03-W10 N11-N14
      4 1262  949 (S01W15 L=276)   583    138        1 (1197/ 916-1319/ 982)  23.2 W11-W18 S03-N01
      5 1703  939 (N02W47 L=308)   647    106        4 (1662/ 902-1754/ 969)  23.2 W43-W51 N00-N03
      6  279  852 (S16E53 L=208)   760     84        2 ( 240/ 813- 313/ 886) -18.4 E57-E49 S18-S13
      7 1706  781 (S08W50 L=311)   293     85        0 (1689/ 750-1723/ 824)  90.0 W49-W51 S10-S05
      8  266  776 (S21E56 L=206)   915    103        1 ( 244/ 729- 287/ 818) -68.2 E58-E52 S24-S18
      9 1232  739 (S14W16 L=277)   986    328        2 (1103/ 671-1359/ 809) -20.2 W06-W25 S17-S10
     10 1017  727 (S17W02 L=263)   310    110        0 ( 971/ 708-1069/ 741)   0.0 W00-W05 S17-S15
     11  466  674 (S26E37 L=224)   565    137        0 ( 406/ 658- 529/ 688)  15.3 E42-E31 S27-S24
  

    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