Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1521 1732 (N34W61 L= 75)  1678    261        4 (1402/1690-1640/1779) -21.4 W50-W74 N29-N39
      2  635 1708 (N56E17 L=357)  1647    151        4 ( 583/1650- 696/1754) -23.2 E24-E10 N50-N59
      3  945 1610 (N39W08 L= 22)  1432    291        1 ( 833/1540-1084/1653) -42.9 W00-W15 N31-N45
      4  369 1618 (N55E53 L=321)   614     80        1 ( 344/1584- 405/1647)  90.0 E57-E46 N51-N56
      5 1077 1409 (N23W11 L= 25)   474     80        1 (1043/1381-1118/1433)   0.0 W08-W14 N21-N25
      6  386 1294 (N30E39 L=335)  1709    382        5 ( 269/1192- 509/1393) -34.2 E54-E29 N21-N38
      7  491 1333 (N30E30 L=344)   274     47        0 ( 461/1313- 514/1351) -11.3 E32-E27 N28-N31
      8 1222 1306 (N14W17 L= 31)   589    142        0 (1148/1267-1278/1348) -49.4 W13-W20 N11-N18
      9  202 1108 (N22E58 L=316)  2736    354        5 (  93/1006- 339/1216) -22.3 E84-E44 N13-N29
     10  429 1099 (N17E36 L=338)   339     49        1 ( 411/1069- 447/1130) -71.6 E37-E35 N15-N18
     11  538 1046 (N12E29 L=345)   282    123        2 ( 515/1015- 578/1095) -63.4 E30-E26 N09-N15
     12  127 1029 (N18E70 L=305)   746     63        1 ( 103/1004- 160/1053) -18.4 E74-E64 N16-N19
     13  481  819 (S00E37 L=337)   479     76        0 ( 440/ 798- 517/ 838) -26.6 E40-E34 S01-N01
     14 1341  644 (S27W12 L= 26)   374     82        2 (1313/ 606-1368/ 683)  63.4 W09-W14 S29-S24
     15  545  533 (S18E42 L=332)   313     35        2 ( 520/ 508- 572/ 552) -45.0 E43-E40 S19-S16
     16  341  466 (S19E66 L=309)  2447    131        5 ( 279/ 419- 387/ 512) -33.7 E73-E59 S21-S15
  

    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