Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  703 1700 (N26E40 L=119)   806    204        1 ( 642/1657- 769/1748) -17.1 E46-E33 N23-N28
      2  461 1593 (N14E52 L=107)   324     72        1 ( 438/1556- 490/1617)  78.7 E54-E50 N11-N16
      3  392 1486 (N07E53 L=106)   854    240        1 ( 326/1422- 478/1538)  28.6 E57-E45 N01-N10
      4  621 1499 (N11E36 L=123)   291    100        0 ( 594/1459- 652/1536) -33.7 E38-E32 N09-N13
      5  911 1455 (N15E17 L=142)   277     87        0 ( 870/1431- 955/1484)  -8.1 E20-E13 N13-N15
      6  585 1364 (N03E34 L=125)   233     74        0 ( 564/1337- 611/1398) -31.0 E36-E31 N01-N04
      7  808  791 (S25E07 L=152)   222     59        0 ( 797/ 756- 821/ 830) -53.1 E08-E05 S27-S23
      8 1744  635 (S09W59 L=218)   394     70        1 (1722/ 603-1766/ 662)  63.4 W56-W60 S10-S06
      9 1591  361 (S27W61 L=220)   561     73        0 (1553/ 342-1623/ 384)  59.0 W58-W63 S29-S24
     10  858  208 (S65W24 L=183)   486     93        3 ( 832/ 181- 888/ 232)  21.8 W19-W28 S68-S62
  

    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