Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1653 (N43E06 L=319)   485     98        1 (1085/1613-1135/1692) -45.0 E08-E03 N39-N45
      2  347 1449 (N15E58 L=267)  1031    147        1 ( 320/1386- 386/1509)  90.0 E60-E55 N10-N19
      3 1671 1340 (N30W43 L=  8)   619    214        1 (1589/1272-1715/1428) -36.0 W35-W47 N26-N34
      4 1019 1268 (N14E04 L=321)   241    104        1 ( 986/1243-1051/1318) -45.0 E06-E01 N12-N17
      5 1785 1242 (N25W54 L= 19)   495     37        1 (1763/1216-1810/1266)  53.1 W51-W57 N23-N27
      6 1185  732 (S15W15 L=340)   854    277        1 (1112/ 645-1235/ 785)  90.0 W09-W17 S20-S11
      7 1286  675 (S17W23 L=347)   269     83        3 (1264/ 641-1311/ 706)  63.4 W20-W23 S19-S15
      8 1392  661 (S16W30 L=355)   418    171        1 (1328/ 605-1447/ 706)  58.0 W27-W32 S20-S12
      9  523  626 (S33E30 L=295)   929    210        2 ( 431/ 588- 616/ 674)  20.6 E38-E22 S37-S28
     10 1556  467 (S26W52 L= 17)  3965    574        6 (1335/ 388-1781/ 616)  35.4 W36-W67 S35-S13
     11  968  384 (S43W09 L=334)  1278    276        3 ( 852/ 351-1093/ 431)  25.2 W00-W18 S47-S37
  

    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