Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  834 1697 (N49E20 L= 56)   346     56        2 ( 804/1672- 862/1716)  -9.5 E23-E17 N46-N50
      2 1559 1565 (N39W46 L=122)   748    134        0 (1506/1529-1614/1609) -36.9 W42-W50 N36-N42
      3  731 1434 (N29E22 L= 54)   288     76        0 ( 716/1403- 752/1471) -59.0 E23-E19 N26-N31
      4  309 1421 (N26E60 L= 16)   658     95        0 ( 278/1378- 339/1464) -31.0 E66-E55 N22-N28
      5 1410 1366 (N25W27 L=103)   291    101        0 (1364/1338-1446/1393)  23.2 W22-W29 N23-N26
      6 1646 1276 (N19W44 L=121)   806    181        1 (1583/1212-1700/1337) -45.0 W39-W48 N15-N23
      7  571 1105 (N07E30 L= 47)   235     75        0 ( 541/1076- 587/1144) -53.1 E31-E28 N05-N09
      8 1752  856 (S07W53 L=129)   566     10        0 (1729/ 832-1775/ 880) -36.9 W50-W55 S08-S05
      9 1293  838 (S08W17 L= 93)   275    118        0 (1246/ 820-1331/ 865)  11.3 W14-W19 S09-S06
     10 1139  721 (S16W08 L= 84)   936    246        1 (1081/ 631-1171/ 831)  77.0 W04-W10 S21-S08
     11  797  734 (S15E14 L= 62)   486    164        0 ( 746/ 694- 868/ 790) -32.0 E17-E09 S17-S11
     12  302  708 (S19E55 L= 22)  1327    190        0 ( 242/ 646- 390/ 747)  14.7 E64-E45 S23-S16
     13 1426  672 (S19W28 L=104)   462     90        0 (1389/ 638-1462/ 711) -33.7 W24-W30 S20-S16
     14 1643  585 (S25W48 L=125)  1513    176        2 (1607/ 505-1689/ 656)  84.3 W44-W51 S30-S19
     15  934  576 (S26E05 L= 71)  1106    393        1 ( 734/ 530-1084/ 617)  14.6 E19-W04 S29-S22
     16 1269  538 (S28W18 L= 94)   371    128        0 (1216/ 507-1317/ 573)  35.5 W14-W21 S30-S25
  

    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