Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1039 1630 (N44W00 L=125)  1170    262        6 ( 981/1555-1095/1697) -85.2 E05-W04 N37-N49
      2 1421 1633 (N44W34 L=159)  1237    210        1 (1381/1550-1460/1696)  50.2 W29-W40 N37-N49
      3  358 1542 (N34E56 L= 69)   498     95        2 ( 337/1514- 382/1566) -36.9 E59-E53 N32-N35
      4  581 1448 (N29E33 L= 92)   269     92        1 ( 565/1405- 593/1488) -68.2 E34-E31 N26-N31
      5 1455 1437 (N30W31 L=156)   334     93        2 (1416/1393-1485/1464) -51.3 W28-W32 N26-N31
      6 1194 1290 (N20W11 L=136)   258     83        0 (1155/1266-1244/1315) -26.6 W08-W14 N18-N21
      7  659 1204 (N14E24 L=101)   383    165        0 ( 632/1137- 692/1280) -66.0 E26-E21 N09-N18
      8 1709 1195 (N14W48 L=173)   289     64        2 (1685/1176-1731/1212)   0.0 W46-W50 N12-N14
      9 1168 1089 (N07W09 L=134)   281    140        1 (1091/1073-1235/1112) -12.5 W04-W13 N06-N08
     10 1482  880 (S05W30 L=155)   290    102        1 (1467/ 831-1499/ 932)  90.0 W28-W30 S08-S02
     11  394  852 (S09E43 L= 82)   312    115        0 ( 354/ 826- 445/ 883) -23.2 E46-E38 S10-S07
     12  352  676 (S21E49 L= 76)   417    108        0 ( 317/ 652- 393/ 713) -26.6 E51-E44 S22-S18
     13  232  625 (S25E63 L= 62)   703    235        2 ( 194/ 608- 262/ 646)  18.4 E68-E59 S27-S23
     14  634  379 (S42E32 L= 93)   997    187        4 ( 591/ 331- 673/ 440)  45.0 E37-E27 S46-S36
     15  538  403 (S40E40 L= 85)   496     80        1 ( 512/ 376- 559/ 439)  51.3 E42-E37 S42-S37
     16 1465  379 (S41W38 L=164)  2036    364        2 (1314/ 338-1610/ 439)  11.3 W25-W51 S43-S35
  

    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