Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1095 1683 (N34E14 L=308)  1177    237        2 (1001/1613-1179/1726) -10.6 E21-E05 N31-N36
      2  350 1548 (N10E60 L=261)   411     82        2 ( 308/1518- 372/1577)  45.0 E63-E57 N07-N11
      3 1346 1405 (N22W10 L=331)   207     67        0 (1329/1372-1362/1441) -78.7 W08-W10 N19-N24
      4 1198 1401 (N18W01 L=322)   392     83        0 (1157/1376-1249/1426)   0.0 E01-W04 N16-N19
      5  311 1389 (S01E56 L=265)   732    121        3 ( 251/1371- 368/1413)   0.0 E63-E50 S01-N00
      6 1865 1103 (N20W56 L= 17)  2744    376        5 (1828/ 998-1936/1262) -86.6 W50-W67 N14-N31
      7 1662 1101 (N13W36 L=358)   669    159        1 (1638/1043-1702/1172) -45.0 W34-W40 N10-N17
      8 1395 1064 (N04W20 L=341)   575    150        1 (1308/1035-1469/1102)  41.2 W15-W23 N00-N07
      9  510  886 (S26E29 L=292)   295     44        0 ( 480/ 863- 539/ 907)   0.0 E31-E26 S27-S24
     10 1806  759 (S00W59 L= 20)   595     57        0 (1768/ 745-1844/ 776)  14.0 W54-W63 S01-N01
     11 1696  557 (S14W59 L= 20)  1674    369        7 (1624/ 446-1748/ 693) -59.0 W52-W64 S21-S06
     12  288  643 (S46E53 L=269)  2534    313        5 ( 160/ 596- 440/ 685)   1.2 E82-E33 S47-S43
     13 1100  628 (S27W15 L=336)   211     76        0 (1063/ 603-1134/ 655) -11.3 W12-W17 S28-S26
     14  565  559 (S46E21 L=300)   257     62        0 ( 540/ 528- 589/ 584) -31.0 E23-E18 S47-S44
     15 1464  480 (S25W45 L=  6)   567     97        1 (1419/ 459-1509/ 504)   0.0 W40-W49 S26-S23
     16 1465  326 (S34W57 L= 19)  2184    285        4 (1322/ 298-1598/ 369)  35.8 W46-W68 S39-S26
  

    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