Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  441 1630 (N19E55 L=276)  1315    297        1 ( 413/1581- 493/1708) -28.1 E62-E47 N16-N25
      2  379 1409 (N05E50 L=280)   526     92        3 ( 351/1374- 402/1444) -53.1 E52-E47 N02-N06
      3 1010 1284 (N08E06 L=324)   560    236        0 ( 952/1216-1105/1319)   0.0 E09-W00 N02-N09
      4  387 1083 (S15E42 L=289)   514    164        1 ( 324/1038- 442/1125)  41.2 E46-E37 S18-S11
      5  936 1038 (S08E05 L=325)   241    112        0 ( 884/1023- 999/1063)   0.0 E08-E01 S08-S07
      6 1089  832 (S17W08 L=338)   183     46        0 (1068/ 807-1106/ 860)  76.0 W06-W08 S18-S14
      7  827  768 (S26E07 L=323)   408    157        0 ( 753/ 729- 902/ 805)  45.0 E10-E02 S29-S21
      8  860  204 (S66W17 L=348)   633     81        2 ( 834/ 175- 884/ 233) -24.4 W11-W23 S69-S63
  

    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