Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  904 1905 (N47E53 L=220)   606     46        1 ( 884/1875- 927/1926) -33.7 E57-E47 N44-N50
      2  780 1832 (N37E51 L=222)   886     94        3 ( 758/1783- 801/1885) -34.7 E58-E45 N32-N41
      3 1077 1636 (N31E17 L=256)   477     98        1 (1052/1580-1090/1687) -58.0 E18-E13 N27-N35
      4 1330 1562 (N35W02 L=275)  1044    290        0 (1261/1432-1389/1666) -40.6 E05-W09 N27-N40
      5 1645 1594 (N52W30 L=303)   408     48        0 (1627/1559-1665/1616)  54.5 W27-W34 N47-N54
      6 1464 1124 (N12W22 L=294)  1370    417        1 (1400/ 955-1505/1289) -73.3 W16-W25 N01-N21
      7  477 1162 (S12E37 L=236)   471     96        1 ( 456/1111- 506/1212) -56.3 E38-E33 S15-S09
      8  632  965 (S20E21 L=251)   213     64        0 ( 607/ 933- 655/ 990) -26.6 E23-E19 S21-S18
      9  733  709 (S32E08 L=264)   215     82        0 ( 712/ 676- 745/ 753)  90.0 E09-E07 S34-S28
     10  984  418 (S42W20 L=292)   626    233        1 ( 879/ 401-1088/ 438)  17.4 W11-W27 S45-S38
     11 1297  411 (S31W40 L=313)   573    126        0 (1242/ 395-1358/ 436)  11.3 W35-W45 S33-S29
  

    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