Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  398 1623 (N22E65 L=197)   684     56        1 ( 375/1595- 427/1654) -14.0 E71-E59 N20-N23
      2  353 1474 (N12E59 L=204)  1159    158        2 ( 328/1416- 394/1547) -66.0 E63-E52 N07-N16
      3  393 1380 (N07E50 L=213)   347    107        0 ( 369/1341- 421/1427)  80.5 E51-E48 N04-N10
      4 1609 1306 (N28W34 L=296)   688    122        1 (1565/1263-1649/1361) -26.6 W29-W37 N25-N30
      5  228  821 (S30E58 L=205)   917    227        0 ( 171/ 748- 289/ 912) -36.9 E65-E50 S33-S24
      6  535  834 (S23E28 L=235)   373    132        0 ( 489/ 793- 570/ 895)  63.4 E30-E25 S26-S18
      7 1032  699 (S21W07 L=270)   425    115        2 ( 982/ 669-1087/ 735)  45.0 W04-W10 S23-S17
      8  688  579 (S36E13 L=249)   651    163        1 ( 621/ 543- 760/ 624) -10.3 E18-E07 S37-S34
      9  816  579 (S33E04 L=259)   216     47        1 ( 798/ 557- 843/ 606)  90.0 E04-E01 S35-S31
     10  279  575 (S47E60 L=202)   977    115        1 ( 223/ 552- 331/ 603)   8.6 E83-E50 S49-S43
     11 1554  480 (S22W52 L=314)   588    113        0 (1526/ 430-1593/ 518) -81.9 W50-W53 S26-S19
     12 1637  435 (S23W64 L=326)  1734    235        5 (1568/ 367-1710/ 518)  80.5 W59-W68 S28-S16
  

    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