Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1108 1873 (N66E14 L=  6)   608     59        1 (1084/1846-1139/1894)  66.8 E16-E10 N62-N70
      2 1631 1595 (N48W50 L= 70)  1932    176        3 (1549/1562-1704/1629)   9.2 W37-W70 N44-N51
      3 1301 1488 (N34W13 L= 33)  1460    398        3 (1234/1350-1390/1565) -52.4 W07-W20 N25-N39
      4  402 1469 (N20E53 L=327)   533     76        1 ( 366/1451- 440/1488)  -6.3 E57-E48 N18-N20
      5  693 1423 (N21E28 L=352)   289     66        2 ( 658/1399- 723/1443)  14.0 E29-E25 N18-N22
      6 1183 1312 (N20W06 L= 27)   351    101        1 (1149/1285-1233/1353) -26.6 W03-W09 N18-N22
      7 1260 1167 (N12W13 L= 33)   241     51        3 (1236/1146-1290/1190) -18.4 W11-W14 N11-N13
      8 1330 1033 (N05W19 L= 39)   330    122        0 (1307/ 980-1355/1090)  74.1 W17-W19 N01-N08
      9 1306  829 (S08W20 L= 40)   251    115        0 (1289/ 781-1321/ 884) -80.5 W18-W21 S10-S04
     10 1816  809 (S03W62 L= 82)   426    111        0 (1803/ 770-1830/ 869) -78.7 W60-W63 S05-N00
     11  333  759 (S25E48 L=333)  1180    113        2 ( 294/ 701- 363/ 823)  69.4 E51-E44 S28-S20
     12 1027  726 (S18W04 L= 24)   245      1        1 (1005/ 702-1059/ 744)  45.0 W02-W05 S19-S15
     13  919  579 (S29E01 L= 19)   812    212        1 ( 831/ 551-1000/ 619)  20.0 E07-W04 S31-S24
     14 1295  570 (S24W25 L= 45)   249    106        0 (1273/ 549-1335/ 606) -21.8 W22-W27 S25-S22
     15  421  567 (S37E42 L=338)   803     98        1 ( 381/ 524- 471/ 602)   0.0 E46-E36 S40-S35
     16 1399  495 (S28W35 L= 55)  3381    760        2 (1075/ 426-1709/ 571)  18.1 W12-W61 S36-S19
  

    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