Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1543 (N26E10 L=197)   851    247        0 ( 831/1472-1019/1594) -19.7 E17-E03 N22-N29
      2  272 1373 (N13E56 L=151)  3560    387        4 ( 196/1212- 369/1543) -51.0 E68-E45 N03-N24
      3  444 1421 (N16E42 L=165)   473     59        2 ( 415/1387- 480/1449)  45.0 E44-E39 N13-N17
      4  358 1408 (N15E49 L=158)   754     80        2 ( 324/1375- 390/1445)  63.4 E52-E45 N13-N17
      5  499 1401 (N15E37 L=170)   491    123        0 ( 447/1361- 535/1445)  51.3 E40-E34 N12-N17
      6 1487 1289 (N15W27 L=234)   384     55        1 (1455/1257-1514/1310) -56.3 W24-W29 N12-N15
      7 1313 1277 (N12W16 L=223)   228     98        0 (1267/1257-1355/1293)  18.4 W12-W18 N10-N13
      8  215 1090 (S04E57 L=150)   546     92        3 ( 177/1060- 258/1124)  20.6 E61-E52 S06-S02
      9  608 1034 (S07E25 L=182)   296     86        3 ( 586/ 991- 625/1073) -78.7 E26-E24 S09-S04
     10  184  968 (S11E61 L=146)   659    111        3 ( 166/ 915- 209/1021)  45.0 E63-E57 S14-S08
     11  521  975 (S10E31 L=176)   261     56        0 ( 501/ 948- 548/ 997)  45.0 E32-E29 S12-S09
     12  362  880 (S17E44 L=164)  1970    407        4 ( 252/ 776- 449/ 948)  32.0 E55-E36 S22-S12
     13  578  907 (S14E27 L=180)   349    146        0 ( 508/ 876- 620/ 947)  26.6 E32-E24 S16-S11
     14  197  855 (S18E60 L=147)  1997    186        5 ( 144/ 786- 264/ 922)  29.7 E67-E52 S22-S14
     15  522  838 (S19E31 L=176)   377    114        1 ( 488/ 805- 577/ 878) -33.7 E33-E27 S20-S16
     16  208  678 (S29E65 L=142)   614     59        0 ( 193/ 641- 223/ 714)  51.3 E68-E62 S31-S26
     17 1324  514 (S33W26 L=233)   429    115        1 (1262/ 491-1374/ 539)  -6.3 W21-W30 S34-S32
     18 1417  240 (S52W54 L=261)  1705    229        1 (1315/ 209-1524/ 268)   3.0 W40-W80 S54-S49
  

    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