Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  878 1742 (N57W16 L= 39)   388     72        0 ( 847/1723- 908/1763) -31.0 W12-W19 N54-N58
      2 1556 1701 (N29W67 L= 90)   908    150        2 (1491/1672-1617/1737) -48.8 W63-W72 N25-N33
      3 1526 1599 (N25W54 L= 77)   953    211        2 (1444/1549-1610/1635) -45.0 W48-W57 N19-N29
      4  534 1566 (N54E26 L=357)   305     42        1 ( 511/1545- 555/1588) -31.0 E28-E23 N51-N55
      5 1358 1540 (N27W37 L= 61)   334    123        0 (1305/1495-1396/1567) -74.1 W35-W38 N23-N30
      6 1357 1435 (N22W32 L= 55)   436    161        1 (1294/1390-1427/1458) -40.6 W28-W35 N17-N24
      7 1138 1261 (N18W13 L= 36)   775    238        2 (1027/1195-1208/1326) -57.5 W07-W16 N12-N23
      8  414 1151 (N28E38 L=345)   439     94        1 ( 369/1121- 455/1176) -20.6 E42-E34 N25-N30
      9  864 1100 (N15E07 L= 16)   272     80        1 ( 843/1051- 879/1135)  53.1 E08-E05 N12-N17
     10 1874  928 (S23W60 L= 83)  2070    312        3 (1835/ 806-1916/1055)  78.7 W53-W68 S31-S16
     11  451  996 (N18E37 L=346)   298     62        0 ( 438/ 959- 463/1035)  78.7 E38-E35 N15-N20
     12 1176  947 (S01W07 L= 30)   849    209        1 (1126/ 854-1226/1032)  41.2 W02-W11 S05-N02
     13 1343  701 (S19W11 L= 34)   431    157        0 (1270/ 660-1390/ 728) -35.5 W06-W13 S23-S17
     14 1245  525 (S28W01 L= 24)   584    194        0 (1162/ 491-1328/ 573)  -5.2 E04-W07 S29-S26
     15 1531  528 (S37W22 L= 45)   639     90        1 (1495/ 486-1561/ 557)  45.0 W18-W24 S40-S35
     16 1466  447 (S42W15 L= 38)   302     55        0 (1433/ 435-1495/ 460) -11.3 W12-W17 S42-S40
     17 1390  405 (S42W08 L= 31)   327     52        1 (1372/ 379-1415/ 439) -90.0 W06-W10 S44-S38
     18 1296  352 (S43E01 L= 22)   302     46        0 (1282/ 331-1318/ 382) -90.0 E02-W00 S45-S40
     19  678  251 (S34E52 L=331)   706     80        1 ( 653/ 221- 710/ 292)  60.3 E55-E49 S37-S30
  

    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