Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1330 1760 (N55W37 L=162)   635     80        2 (1281/1727-1366/1785)   8.7 W31-W45 N51-N56
      2 1299 1556 (N38W23 L=148)   298     44        1 (1276/1535-1331/1573) -14.0 W21-W25 N36-N39
      3  547 1429 (N30E35 L= 89)  1965    420        2 ( 406/1289- 679/1534) -27.3 E53-E22 N20-N36
      4 1038 1258 (N18W01 L=126)   467    211        0 ( 995/1180-1116/1310)  45.0 E01-W06 N13-N21
      5  232 1225 (N16E61 L= 63)  1055    155        2 ( 211/1152- 256/1301) -60.9 E64-E57 N11-N20
      6  724  976 (N01E19 L=106)  1468    408        3 ( 640/ 828- 849/1104)  47.0 E24-E10 S07-N09
      7  601  638 (S20E30 L= 95)  1385    312        0 ( 516/ 507- 680/ 733)  43.2 E39-E22 S29-S13
      8 1479  571 (S27W32 L=157)  2409    412        4 (1338/ 449-1583/ 730)  52.1 W22-W38 S35-S16
      9  919  668 (S18E07 L=117)   235     49        1 ( 901/ 647- 940/ 693)  45.0 E08-E06 S19-S16
     10 1698  643 (S23W51 L=175)   342     67        0 (1685/ 606-1710/ 675) -78.7 W49-W52 S25-S20
     11  957  476 (S32E06 L=119)  1578    505        0 ( 813/ 369-1180/ 633)  35.1 E18-W09 S40-S21
     12 1197  404 (S38W12 L=137)   849     96        1 (1156/ 359-1246/ 443)  26.6 W09-W16 S41-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