Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  603 1808 (N31E56 L=123)   682    128        1 ( 576/1785- 634/1836) -24.0 E60-E51 N29-N33
      2 1273 1681 (N40E04 L=175)   536     77        0 (1240/1647-1305/1714) -23.2 E07-W00 N38-N42
      3  640 1647 (N20E42 L=137)   350     63        0 ( 608/1633- 673/1662)  11.3 E44-E39 N18-N21
      4  712 1527 (N14E32 L=147)   518    142        0 ( 684/1470- 743/1594) -45.0 E36-E28 N10-N17
      5  371 1008 (S23E41 L=138)   370    122        1 ( 331/ 974- 422/1031)  29.7 E44-E37 S25-S21
      6 1459 1007 (N03W25 L=204)   606    181        2 (1381/ 983-1551/1031)  20.0 W20-W31 N00-N05
      7 1572  889 (N01W35 L=215)   293     75        1 (1541/ 869-1595/ 920)  45.0 W33-W37 S01-N02
      8  709  834 (S26E14 L=165)   516    187        0 ( 614/ 803- 789/ 867)  32.5 E20-E09 S29-S21
      9 1597  773 (S05W41 L=220)   664    102        1 (1568/ 724-1617/ 834) -63.4 W38-W43 S07-S01
     10 1340  693 (S17W28 L=207)  1644    350        0 (1171/ 639-1504/ 775)  43.2 W18-W35 S22-S07
     11 1539  708 (S10W40 L=219)   299    100        0 (1515/ 659-1553/ 743) -38.7 W36-W41 S12-S08
     12 1015  455 (S40W18 L=197)  1797    361        5 ( 857/ 424-1186/ 488)  16.3 W06-W30 S44-S35
     13 1520  395 (S27W54 L=233)  1688    380        1 (1359/ 364-1652/ 425)  31.0 W43-W64 S32-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