Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1312 1829 (N44W56 L=206)   967    115        2 (1270/1777-1352/1866)  13.2 W48-W65 N40-N45
      2 1510 1556 (N21W47 L=197)  1383    219        3 (1473/1439-1538/1656)  32.9 W38-W55 N15-N26
      3  928 1526 (N35W07 L=157)   872    237        1 ( 850/1443-1027/1580)  14.9 W00-W15 N31-N37
      4 1285 1563 (N28W32 L=181)   300    109        0 (1239/1538-1332/1580)  -8.1 W27-W34 N25-N28
      5 1361  944 (S09W17 L=167)   208     45        1 (1337/ 922-1380/ 974)  33.7 W15-W18 S09-S07
      6 1400  929 (S11W19 L=169)   285     63        0 (1376/ 892-1424/ 957)  36.9 W17-W21 S12-S09
      7 1689  717 (S31W38 L=188)   408     72        0 (1668/ 672-1708/ 756)  59.0 W35-W39 S34-S28
      8 1159  594 (S25E02 L=147)   361    136        0 (1112/ 535-1208/ 651)  35.5 E06-W01 S27-S22
      9  855  407 (S29E28 L=122)  1290    316        1 ( 820/ 253- 881/ 550)  56.3 E37-E23 S40-S19
     10  568  436 (S21E46 L=104)   565    143        0 ( 526/ 384- 608/ 489)  22.6 E52-E40 S23-S17
     11  439  438 (S18E57 L= 93)   402    101        0 ( 405/ 415- 458/ 485)  20.6 E61-E52 S19-S15
     12  820  236 (S41E42 L=108)   387     33        1 ( 798/ 217- 846/ 257)   9.5 E45-E39 S41-S38
  

    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