Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1077 1784 (N42E25 L=201)   633    184        0 (1029/1699-1122/1849) -29.1 E34-E16 N36-N46
      2  569 1697 (N21E54 L=172)   378     45        0 ( 550/1676- 593/1714)  63.4 E56-E51 N19-N22
      3 1019 1440 (N17E12 L=213)   775    207        0 ( 950/1349-1089/1501) -22.6 E17-E05 N13-N18
      4  236 1464 (N00E73 L=153)   820     96        1 ( 217/1428- 264/1495)  90.0 E78-E67 S02-N02
      5 1867 1045 (N22W58 L=283)  1170    186        1 (1807/ 959-1896/1097) -22.6 W49-W63 N17-N24
      6  617  990 (S19E23 L=203)   860    366        0 ( 479/ 896- 769/1042)  43.2 E30-E14 S27-S12
      7  281  838 (S37E49 L=176)   372     80        0 ( 266/ 798- 305/ 875)  56.3 E51-E46 S39-S33
      8  744  831 (S25E11 L=215)   525    167        0 ( 682/ 787- 824/ 861)  48.8 E14-E06 S29-S21
      9  908  661 (S30W04 L=230)  1804    569        3 ( 679/ 563-1149/ 716)  12.0 E11-W22 S35-S26
     10 1452  526 (S20W43 L=269)   592    140        1 (1399/ 486-1519/ 569)  63.4 W40-W45 S23-S15
     11  346  478 (S62E53 L=173)   537     43        1 ( 320/ 460- 370/ 504)  10.6 E62-E46 S63-S59
     12 1534  322 (S27W67 L=293)   635    115        1 (1485/ 287-1578/ 357)  90.0 W64-W68 S31-S24
  

    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