Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  550 1600 (N15E48 L=262)   359     85        1 ( 531/1569- 578/1627) -36.9 E50-E45 N13-N16
      2  331 1536 (N07E64 L=246)  2144    206        2 ( 262/1464- 379/1617)  84.8 E71-E61 N01-N12
      3  786 1564 (N19E30 L=280)   290    102        1 ( 748/1532- 825/1600) -14.0 E34-E26 N17-N20
      4  968 1410 (N14E14 L=296)   372    109        1 ( 935/1363- 995/1464) -81.9 E15-E12 N11-N18
      5 1545 1068 (N12W29 L=339)   405    105        0 (1519/1023-1581/1118) -23.2 W25-W32 N10-N13
      6 1659 1084 (N16W36 L=347)   384    102        1 (1629/1047-1698/1111)  51.3 W35-W39 N13-N18
      7 1705  873 (N06W46 L=356)  1947    337        5 (1627/ 780-1748/ 960) -21.8 W37-W52 N01-N10
      8 1852  825 (N09W62 L= 12)   630    134        3 (1824/ 784-1875/ 857) -14.0 W57-W66 N07-N10
      9  553  695 (S37E21 L=289)   570    146        1 ( 498/ 659- 627/ 726)  37.9 E25-E16 S40-S33
     10 1686  439 (S16W67 L= 17)  1093    243        1 (1666/ 383-1698/ 504) -24.4 W62-W73 S18-S12
  

    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