Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  761 1705 (N34E32 L=255)   935    189        0 ( 714/1611- 786/1781) -47.3 E40-E27 N26-N39
      2  467 1629 (N25E52 L=234)   411     48        1 ( 441/1612- 495/1650)  26.6 E54-E49 N23-N26
      3 1214 1612 (N34W03 L=289)   255     71        0 (1188/1574-1235/1643) -53.1 W01-W05 N31-N35
      4  356 1537 (N18E58 L=228)   887    141        2 ( 306/1488- 378/1573) -63.4 E64-E54 N14-N20
      5 1915 1107 (N15W67 L=353)   680     75        1 (1893/1065-1933/1141) -53.1 W63-W70 N13-N17
      6  161  993 (S17E65 L=221)  2382    313        2 ( 110/ 875- 236/1134)  59.0 E75-E56 S23-S08
      7  881  864 (S18E06 L=280)   314    137        0 ( 812/ 849- 946/ 885)   7.1 E10-E02 S18-S16
      8  498  812 (S26E32 L=254)   996    233        3 ( 402/ 773- 617/ 859)  16.4 E40-E23 S28-S21
      9  830  813 (S21E09 L=277)   316    126        0 ( 784/ 773- 887/ 852) -20.6 E12-E04 S22-S19
     10 1520  601 (S23W41 L=327)  1493    433        1 (1355/ 490-1659/ 745)  58.2 W33-W47 S32-S11
     11  330  710 (S34E49 L=237)   422     43        0 ( 310/ 684- 348/ 740)  76.0 E51-E47 S35-S31
     12  725  289 (S61E15 L=272)   826    148        2 ( 661/ 256- 778/ 318)  -3.6 E23-E07 S63-S58
  

    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