Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  660 1592 (N47E24 L=280)   317     70        3 ( 621/1576- 700/1610)   0.0 E27-E19 N45-N48
      2 1161 1481 (N32W16 L=319)  1209    305        3 (1025/1428-1300/1523) -27.9 W06-W24 N26-N36
      3 1695 1470 (N23W59 L=  2)   470     85        1 (1651/1447-1733/1494) -35.5 W55-W62 N20-N25
      4  606 1288 (N26E26 L=278)   261     95        1 ( 564/1275- 657/1298)  -8.1 E28-E21 N25-N27
      5  716 1172 (N18E18 L=285)   302    100        1 ( 668/1146- 751/1207)   9.5 E21-E15 N16-N20
      6  387 1115 (N17E43 L=261)   291     50        2 ( 360/1095- 419/1129)  11.3 E45-E40 N15-N17
      7 1496 1004 (S01W29 L=333)   266    135        1 (1450/ 963-1532/1061) -63.4 W26-W31 S04-N02
      8 1799  835 (S17W54 L=358)  5821    603       14 (1679/ 611-1935/1015)  42.5 W43-W75 S30-S08
      9 1378  841 (S10W20 L=323)   715    164        4 (1330/ 784-1417/ 891)  56.3 W16-W22 S13-S06
     10  179  760 (S05E71 L=233)  3501    267        8 ( 139/ 673- 217/ 875) -85.2 E87-E66 S10-N02
     11 1533  842 (S12W31 L=334)   388    118        1 (1491/ 799-1562/ 872) -51.3 W27-W32 S14-S09
     12 1665  776 (S18W41 L=344)   387    113        1 (1632/ 728-1705/ 806)  33.7 W38-W44 S21-S17
     13 1464  714 (S19W25 L=328)   570    195        3 (1386/ 664-1538/ 757)  15.3 W19-W30 S21-S16
     14  844  694 (S13E15 L=288)   386     97        0 ( 796/ 673- 891/ 717) -21.8 E17-E12 S14-S10
     15  802  533 (S23E21 L=282)  1505    377        3 ( 739/ 443- 923/ 676)  86.4 E24-E13 S29-S13
     16  400  534 (S20E53 L=250)   796    149        3 ( 355/ 503- 449/ 573) -23.2 E57-E48 S22-S16
     17  611  465 (S26E37 L=266)   384     75        1 ( 592/ 430- 636/ 497)  38.7 E39-E34 S28-S23
     18 1594  406 (S45W43 L=346)  1546    186        2 (1544/ 351-1626/ 491) -51.8 W38-W50 S51-S37
     19 1045  312 (S43E09 L=294)   570     54        0 (1015/ 285-1079/ 334)  31.0 E11-E06 S45-S41
  

    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