Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  822 1605 (N36E12 L=184)   451     81        2 ( 775/1586- 868/1626)   0.0 E15-E08 N34-N37
      2 1598 1563 (N30W45 L=241)   509    115        1 (1567/1530-1626/1600)  23.2 W41-W48 N27-N32
      3  755 1303 (N16E16 L=180)   303    108        0 ( 709/1282- 818/1324) -15.9 E18-E11 N14-N17
      4  149 1235 (N15E66 L=130)   923    183        0 ( 115/1202- 182/1271) -33.7 E71-E60 N12-N18
      5 1395 1183 (N06W23 L=219)   202     32        2 (1376/1165-1413/1205)  56.3 W22-W24 N04-N07
      6 1050 1011 (S03W01 L=197)   638    155        2 (1024/ 938-1070/1108) -90.0 W00-W02 S07-N02
      7  219  956 (S02E57 L=139)   329     57        0 ( 194/ 936- 241/ 973) -26.6 E59-E54 S03-N00
      8 1093  838 (S14W03 L=199)   265     51        0 (1076/ 806-1113/ 870)  63.4 W02-W04 S15-S11
      9  684  514 (S32E27 L=169)   488     70        1 ( 651/ 489- 723/ 551) -45.0 E29-E24 S34-S29
     10  559  516 (S32E37 L=159)   527    115        1 ( 516/ 485- 612/ 538)  14.0 E40-E32 S33-S29
     11  502  448 (S36E44 L=152)   496     94        0 ( 470/ 427- 542/ 472)   8.1 E47-E40 S37-S34
  

    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