Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  349 1350 (N35E49 L= 65)   368     47        0 ( 321/1333- 380/1365)   0.0 E52-E45 N33-N35
      2 1188 1257 (N18W14 L=128)   380    166        0 (1166/1193-1230/1326) -83.7 W12-W15 N13-N22
      3 1275 1285 (N18W20 L=134)   378    137        1 (1218/1231-1316/1325) -60.3 W17-W22 N14-N21
      4  868 1111 (N14E08 L=106)   214     68        1 ( 840/1090- 900/1141)  14.0 E10-E06 N12-N15
      5 1692 1028 (S05W44 L=158)   514    117        0 (1644/1002-1738/1060) -45.0 W40-W48 S07-S01
      6  951  625 (S17E10 L=104)   565    187        3 ( 897/ 555- 998/ 703)  45.0 E15-E06 S21-S12
      7 1195  524 (S28W04 L=118)   795    319        1 (1097/ 418-1291/ 634)  42.9 E04-W11 S34-S21
      8  609  481 (S23E38 L= 76)   772    239        1 ( 524/ 401- 681/ 546)  22.8 E49-E30 S27-S18
      9  739  403 (S30E31 L= 83)  1006    258        2 ( 642/ 343- 838/ 481)  20.9 E42-E21 S33-S25
     10 1199  396 (S38W02 L=116)   492    158        1 (1126/ 364-1260/ 438)  21.8 E03-W07 S39-S35
  

    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