Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1054 1722 (N53W08 L= 53)   363     83        0 (1013/1706-1106/1732)   0.0 W04-W13 N51-N53
      2 1385 1545 (N36W32 L= 77)  1609    543        0 (1193/1473-1615/1631) -25.8 W18-W49 N28-N43
      3  391 1425 (N32E49 L=355)   855    154        3 ( 368/1372- 429/1505) -36.9 E55-E43 N28-N37
      4 1147 1392 (N27W10 L= 55)   379    109        2 (1107/1358-1191/1409)  18.4 W07-W13 N25-N28
      5  806 1091 (N10E13 L= 31)   530    113        2 ( 760/1042- 841/1135) -38.7 E16-E11 N06-N12
      6  516 1084 (N10E33 L= 11)   278     93        0 ( 489/1047- 547/1125) -45.0 E35-E30 N07-N12
      7 1826  853 (S12W60 L=104)  2734    493        3 (1705/ 670-1906/1036)  47.9 W49-W71 S22-S01
      8  254  745 (S11E60 L=345)   473     64        0 ( 222/ 719- 286/ 766)  16.7 E65-E55 S13-S09
      9  240  633 (S19E67 L=337)   717     84        1 ( 230/ 589- 257/ 679)  54.5 E71-E64 S22-S15
     10 1465  619 (S24W29 L= 73)   295     67        0 (1434/ 593-1498/ 643)  21.8 W26-W31 S25-S22
     11 1263  551 (S27W14 L= 59)   386    104        0 (1226/ 530-1306/ 575)  26.6 W11-W17 S28-S25
  

    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