Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  736 1758 (N29E45 L=214)  1501    213        4 ( 703/1660- 769/1843) -43.0 E55-E38 N22-N36
      2  928 1710 (N31E28 L=230)   684    150        1 ( 892/1643- 958/1763) -35.0 E33-E23 N26-N34
      3  452 1679 (N18E63 L=196)   763    106        2 ( 414/1640- 492/1713) -71.6 E66-E58 N14-N20
      4 1346 1543 (N33W04 L=263)   605    196        0 (1327/1442-1366/1623) -77.9 W01-W07 N25-N39
      5  342 1590 (N10E67 L=191)  1141     84        3 ( 308/1552- 373/1622) -59.0 E74-E62 N07-N12
      6  715 1584 (N17E35 L=223)   575    124        2 ( 668/1538- 756/1622)  78.7 E37-E32 N13-N20
      7 1737 1402 (N38W40 L=298)  1042    129        1 (1702/1344-1763/1449)  60.9 W35-W43 N33-N42
      8  574 1331 (S01E35 L=224)   366     58        1 ( 551/1295- 599/1359) -26.6 E36-E32 S02-N00
      9  147 1229 (S14E73 L=185)  1868    211        2 ( 114/1118- 171/1308)  90.0 E78-E68 S21-S09
     10  210 1217 (S15E63 L=196)   493     66        3 ( 186/1178- 227/1247)  59.0 E65-E60 S17-S12
     11  369 1171 (S14E45 L=213)   267     63        1 ( 339/1146- 394/1192)   0.0 E48-E42 S15-S13
     12 1679 1089 (N16W38 L=297)   672    208        1 (1600/1053-1752/1141)  30.3 W31-W43 N12-N21
     13  330 1024 (S24E46 L=213)   346     87        0 ( 291/1005- 380/1038)   7.1 E49-E41 S24-S22
     14  264  990 (S27E52 L=206)   339     75        2 ( 227/ 968- 293/1011)  -7.1 E56-E48 S28-S26
     15 1699  888 (N06W45 L=304)   405    160        1 (1685/ 812-1718/ 966) -58.0 W41-W47 N01-N09
     16  784  742 (S29E06 L=252)   229     67        1 ( 757/ 728- 828/ 761)  11.3 E08-E03 S30-S28
     17 1663  611 (S10W54 L=313)   429    130        0 (1596/ 591-1713/ 627)  21.8 W48-W58 S11-S07
     18 1454  439 (S26W49 L=307)  1898    405        3 (1287/ 416-1628/ 471)  19.0 W35-W64 S31-S18
     19 1159  397 (S38W32 L=291)   330     60        1 (1124/ 372-1192/ 416)   9.5 W28-W34 S39-S37
  

    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