Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1050 1877 (N57W49 L=315)   787     79        2 (1008/1845-1082/1907) -68.2 W44-W56 N53-N59
      2  762 1701 (N55W05 L=272)   297     54        0 ( 730/1686- 789/1719)  -8.1 W01-W08 N53-N56
      3 1716 1580 (N15W70 L=337)   562     64        1 (1680/1548-1733/1607)  31.0 W65-W74 N12-N16
      4 1597  782 (S24W29 L=295)   682    161        1 (1516/ 734-1666/ 827)   0.0 W21-W34 S25-S23
      5 1302  684 (S21W07 L=274)   334    153        0 (1277/ 627-1354/ 719)  38.7 W04-W10 S24-S19
      6  426  679 (N00E46 L=220)   277     54        0 ( 409/ 648- 440/ 707)  71.6 E47-E45 S01-N02
      7 1366  620 (S27W09 L=276)   307     93        0 (1344/ 574-1386/ 669)  51.3 W06-W11 S29-S24
      8 1440  561 (S33W13 L=280)   433    100        0 (1388/ 545-1485/ 589) -15.9 W09-W16 S34-S31
      9 1703  508 (S49W39 L=306)  1005    135        1 (1666/ 468-1729/ 568) -50.7 W34-W45 S53-S42
     10 1087  494 (S27E11 L=255)   494    131        2 (1058/ 439-1121/ 550)  40.6 E14-E07 S29-S23
     11  500  349 (S21E60 L=207)   600     50        2 ( 473/ 323- 528/ 378)  53.1 E62-E57 S23-S19
  

    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