Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1305 1752 (N56W14 L= 11)   336     39        0 (1279/1734-1330/1765)  26.6 W11-W17 N53-N58
      2 1111 1572 (N37E02 L=355)   446    190        0 (1050/1511-1185/1626) -32.5 E07-W04 N33-N40
      3  361 1339 (N11E51 L=306)   908    138        3 ( 342/1263- 383/1408) -66.0 E54-E47 N06-N15
      4  705 1287 (N12E24 L=333)   269    139        0 ( 685/1227- 721/1358) -58.0 E26-E21 N08-N16
      5 1779 1333 (N30W56 L= 53)   434     64        1 (1748/1311-1809/1357) -18.4 W52-W60 N28-N30
      6  943 1092 (N04E06 L=351)  1024    386        0 ( 876/ 927-1011/1258)  84.3 E08-E03 S06-N14
      7 1668 1079 (N12W43 L= 40)   358     51        2 (1638/1063-1700/1100) -11.3 W40-W45 N10-N13
      8 1207  630 (S21W17 L= 14)  2264    634        3 ( 960/ 517-1458/ 753)  34.6 W02-W31 S32-S10
      9  469  550 (S38E37 L=319)   634    185        0 ( 422/ 482- 526/ 630) -59.0 E40-E32 S42-S32
     10  646  588 (S32E21 L=336)   337    109        1 ( 588/ 572- 690/ 602)   7.1 E25-E17 S32-S30
  

    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