Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  996 1910 (N63W51 L=199)   533     47        0 ( 983/1882-1009/1943)  13.4 W43-W64 N60-N65
      2  437 1464 (N47E33 L=116)   407     49        1 ( 416/1437- 463/1486)  33.7 E35-E29 N44-N48
      3 1033 1334 (N22W09 L=157)   212     87        0 ( 991/1310-1073/1348)  -9.5 W05-W11 N20-N22
      4  862 1297 (N24E02 L=146)   214     57        1 ( 848/1267- 882/1324) -90.0 E03-W00 N22-N25
      5  340 1117 (N26E43 L=105)   613    175        0 ( 310/1046- 404/1166)  35.5 E46-E37 N22-N28
      6 1698  911 (S20W39 L=188)   473    129        1 (1661/ 864-1753/ 951)  12.5 W35-W44 S21-S18
      7 1623  682 (S32W30 L=178)   524    187        0 (1541/ 635-1692/ 732)   8.1 W22-W36 S32-S30
      8 1254  600 (S26W02 L=151)   207     52        0 (1228/ 565-1277/ 626)  26.6 W00-W04 S27-S24
      9  855  571 (S17E22 L=127)   247     87        0 ( 812/ 543- 890/ 604)  15.9 E25-E18 S17-S15
     10  445  522 (S10E52 L= 97)   815    147        0 ( 421/ 450- 472/ 586)  48.8 E55-E46 S14-S06
     11 1424  499 (S38W11 L=160)   627    146        0 (1356/ 459-1481/ 545)  10.3 W05-W16 S38-S35
     12  897  308 (S35E32 L=117)   341     73        0 ( 872/ 278- 916/ 342)  29.7 E35-E28 S36-S32
     13  559  284 (S28E60 L= 88)  1689    137        3 ( 502/ 241- 618/ 337) -66.8 E64-E54 S31-S23
  

    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