Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1027 1904 (N52E43 L= 90)  1270    148        1 ( 974/1855-1060/1942) -24.2 E55-E35 N46-N56
      2 1093 1830 (N47E28 L=105)   521     75        1 (1079/1795-1112/1869) -37.9 E33-E24 N44-N51
      3  693 1783 (N30E52 L= 82)   468    108        0 ( 661/1735- 734/1817) -19.7 E58-E44 N27-N32
      4 1375 1642 (N43W03 L=136)   302     63        1 (1342/1619-1401/1669)   0.0 W00-W05 N40-N44
      5  710 1625 (N20E38 L= 95)   302     42        1 ( 681/1605- 736/1642)  45.0 E40-E37 N18-N21
      6  844 1496 (N16E24 L=109)   264     85        0 ( 796/1476- 875/1519)  53.1 E26-E22 N13-N17
      7 1447  903 (S00W27 L=161)  1199    373        1 (1347/ 765-1506/1065) -32.9 W16-W34 S06-N05
      8 1500  357 (S27W60 L=193)   741    166        0 (1427/ 330-1578/ 382)  41.6 W55-W65 S30-S22
      9 1292  341 (S35W46 L=179)   721    107        0 (1235/ 312-1330/ 373)  71.6 W43-W48 S39-S32
     10  772  268 (S63W16 L=149)   414     64        0 ( 756/ 241- 799/ 300) -21.8 W11-W21 S64-S59
     11  883  201 (S63W35 L=169)   486     46        2 ( 854/ 179- 909/ 227)  36.9 W31-W40 S64-S59
  

    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