Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1052 1830 (N58W36 L=129)   536    140        4 (1032/1790-1075/1861)  26.6 W31-W41 N55-N60
      2 1146 1665 (N43W29 L=121)   380     67        1 (1122/1636-1180/1691)   9.5 W25-W31 N41-N44
      3 1508 1442 (N19W42 L=135)   968    201        1 (1428/1402-1610/1479) -42.0 W36-W48 N14-N23
      4 1771 1225 (S00W55 L=148)  2274    332        2 (1678/1079-1818/1351) -90.0 W48-W59 S08-N09
      5 1100 1208 (N16W09 L=101)   387    121        0 (1050/1162-1134/1263)  23.2 W04-W12 N13-N18
      6  188 1100 (N26E60 L= 33)   594    148        0 ( 157/1070- 224/1136) -38.7 E63-E56 N24-N28
      7 1804  951 (S17W51 L=144)  1204    164        3 (1751/ 906-1841/1002)  -5.7 W45-W55 S19-S14
      8  705  582 (S13E30 L= 63)  1241    344        0 ( 653/ 452- 786/ 718)  43.2 E37-E20 S20-S05
      9 1545  518 (S38W24 L=117)  1850    272        1 (1487/ 458-1615/ 581) -45.0 W20-W31 S42-S31
     10  516  348 (S25E55 L= 38)   859    254        2 ( 489/ 277- 543/ 417)  61.4 E60-E51 S31-S20
  

    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