Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  618 1529 (N33E34 L=332)   484    100        0 ( 600/1472- 637/1580) -49.4 E37-E30 N29-N36
      2  899 1479 (N31E11 L=354)   489    197        0 ( 816/1419- 972/1511) -22.6 E17-E05 N27-N33
      3 1609 1476 (N33W45 L= 50)   754    118        1 (1548/1440-1655/1507) -24.0 W40-W49 N30-N35
      4 1348 1430 (N29W22 L= 27)   350    186        0 (1262/1377-1407/1461) -21.8 W15-W25 N26-N31
      5  475 1303 (N17E39 L=326)   310     86        1 ( 444/1272- 514/1334) -23.2 E42-E35 N15-N19
      6  153 1266 (N12E76 L=290)  1151    104        1 ( 132/1212- 169/1312) -63.4 E81-E72 N09-N15
      7 1128 1193 (N13W06 L= 11)   250     91        0 (1090/1179-1176/1216) -11.3 W03-W08 N12-N14
      8  477  864 (S10E36 L=330)   912    294        4 ( 428/ 732- 522/ 960)  78.7 E39-E32 S18-S03
      9 1705  905 (S03W47 L= 52)   841    122        3 (1657/ 865-1755/ 951) -24.0 W42-W52 S05-N00
     10 1081  799 (S11W04 L= 10)   261    106        0 (1039/ 762-1115/ 835) -45.0 W01-W06 S13-S08
     11 1434  769 (S12W28 L= 33)   218     66        0 (1404/ 752-1466/ 782)  14.0 W25-W29 S12-S10
     12  997  693 (S18E01 L=  5)   602    122        4 ( 951/ 641-1043/ 738) -45.0 E03-W02 S21-S15
     13  884  653 (S21E08 L=357)   682    155        3 ( 828/ 592- 935/ 702) -35.5 E11-E04 S25-S18
     14  431  560 (S30E44 L=321)  1521    242        1 ( 381/ 500- 526/ 650) -68.2 E47-E35 S34-S24
     15 1281  536 (S28W20 L= 25)  2107    393        4 (1153/ 421-1441/ 645)  35.8 W12-W31 S36-S20
     16  779  581 (S27E15 L=350)   642    138        1 ( 725/ 545- 850/ 613)  29.1 E19-E10 S29-S23
     17 1519  494 (S30W40 L= 46)  1395    273        4 (1401/ 431-1624/ 553)  28.1 W31-W48 S35-S25
     18 1006  413 (S38W01 L=  7)  1005    120        3 ( 955/ 368-1049/ 461)  60.3 E02-W04 S41-S34
  

    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