Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1553 1561 (N35W48 L= 79)  1653    334        1 (1398/1523-1684/1632) -26.6 W35-W62 N30-N41
      2  532 1427 (N32E35 L=355)   895    193        3 ( 483/1367- 598/1508) -34.7 E42-E28 N28-N37
      3 1375 1268 (N18W25 L= 55)   353    164        0 (1321/1228-1443/1316) -40.6 W21-W28 N14-N20
      4  525 1200 (N17E33 L=358)   380    140        0 ( 492/1144- 579/1246)  54.5 E35-E28 N13-N20
      5  697 1120 (N12E20 L= 10)   570    248        1 ( 611/1034- 762/1203) -47.7 E26-E16 N06-N17
      6 1030 1100 (N09W01 L= 31)   468    123        0 ( 994/1045-1059/1143) -90.0 E01-W02 N05-N12
      7 1671 1094 (N04W44 L= 75)   269     58        0 (1642/1077-1701/1111)   9.5 W41-W47 N03-N05
      8 1378  783 (S12W21 L= 52)   452    117        1 (1342/ 748-1438/ 822)  23.2 W18-W25 S14-S10
      9  384  721 (S12E47 L=344)  1077    196        1 ( 351/ 633- 447/ 776)  42.3 E51-E40 S18-S08
  

    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