Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1666 (N25E58 L=148)   463     59        0 ( 424/1639- 464/1691) -76.0 E60-E56 N23-N27
      2 1595 1565 (N40W33 L=239)   338     40        0 (1569/1546-1622/1583)  11.3 W30-W35 N38-N40
      3  330 1554 (N16E63 L=144)   502     51        2 ( 303/1527- 356/1579)  18.4 E65-E59 N14-N18
      4 1423 1427 (N25W18 L=224)   293     85        0 (1384/1401-1462/1454)  -9.5 W15-W21 N24-N26
      5 1050 1435 (N18E06 L=201)   209     30        0 (1025/1418-1075/1451) -14.0 E07-E03 N17-N19
      6 1622 1341 (N24W34 L=240)   424    123        0 (1576/1292-1664/1382) -29.7 W30-W37 N21-N25
      7  136 1110 (S12E70 L=136)  2840    182        5 (  94/1045- 202/1195)  33.7 E79-E61 S16-S06
      8  322  807 (S29E48 L=159)  1321    205        2 ( 229/ 772- 424/ 835)   8.1 E58-E37 S32-S27
      9 1831  674 (S09W68 L=274)   612     73        1 (1793/ 648-1865/ 695)  29.7 W64-W72 S11-S07
     10 1716  644 (S14W56 L=262)   345     38        0 (1694/ 621-1737/ 664)  11.3 W53-W58 S15-S12
     11 1468  509 (S28W42 L=248)  2771    510        4 (1258/ 451-1713/ 580)  28.7 W28-W59 S35-S17
  

    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